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Abstract 

A microscopic approach is developed to account for the magnitudes of the su- 
percurrents observed experimentally in carbon nanotubes placed between supercon- 
ducting contacts. We build up a model for the nanotube ropes encompassing the 
electron repulsion from the Coulomb interaction and the effective attraction given by 
phonon exchange. We show that the available experimental data are consistent with 
the expected decay of the supercurrents along the length of the nanotube samples. 
Our results stress that the propagation of the Cooper pairs is favored in the thick 
ropes, as a consequence of the reduction in the strength of the Coulomb interaction 
from the electrostatic coupling between the metallic nanotubes. We also provide an 
explanation for the temperature dependence of the supercurrents observed in the 
experiments, remarking the existence of a crossover from a very flat behavior at low 
temperatures to a pronounced decay in the long nanotubes. 
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1 Introduction 



Since their discovery in 1991, the carbon nanotubes have shown remarkable electronic 
properties. They may display quite different behaviors, depending on their geometry and 
the quality of the contacts used to probe the electronic transport. It was first predicted[l, 
13 El, and later checked experimentally "U Ej, that the nanotubes have semiconducting or 
metallic behavior depending on the helicity of the hexagonal carbon rings wrapped around 
the tubule. This is the kind of versatility that points at the carbon nanotubes as ideal 
molecules to build electronic devices at the nanometer scale jH]- 

The electron correlations are moreover very strong in the carbon nanotubes, what leads 
to the breakdown of the conventional picture of electronic transport at low energies. The 
Coulomb interaction prevails in the individual nanotubesjTHH], and it drives the electron 
system to a state with the properties of the so-called Luttinger liquid, characterized by 
the absence of quasiparticles at the Fermi level. This feature reflects in the power-law 
behavior of observables like the tunneling density of states, whose suppression at low 
energies has been actually observed in the measurements of the conductance in the carbon 
nanotubes [9*| HU]. 

There have been also experiments revealing the existence of superconducting corre- 
lations in the nanotubes |11| IT"") . The first observations were made in nanotube ropes 
suspended between superconducting electrodes, the most remarkable signal being the ap- 
pearance of supercurrents for temperatures below the critical value of the contacts |llj. 
What has been measured in that experiment is the proximity effect, by which Cooper pairs 
are formed in the nanotubes near the superconducting contacts. Later, superconducting 
transitions have been observed in ropes suspended between metallic, nonsuperconduct- 
ing electrodes [HJ]. The measurements reported in Ref. 03 point at the existence of a 
superconducting phase intrinsic to the carbon nanotubes. More recently, strong supercon- 
ducting correlations have been also reported in individual nanotubes of very short radius, 
inserted in a zeolite matrix [TJ]. 

The observation of superconducting correlations in the nanotubes implies that there is 
a regime which is different to that marked by the dominance of the Coulomb interaction. 
The appearance of the superconducting effects requires a suitable attenuation of the repul- 
sive interaction within the nanotubes, due to either the electrostatic coupling among the 
metallic nanotubes in a rope or the coupling with nearby conductors. Another factor that 
seems to be crucial is the achievement of high quality contacts between the nanotubes and 
the electrodes. This has been possible in the experiments described in Ref. [H] by the 
use of a technique that allows to sold and suspend the nanotubes between the contacts. 
The electrodes were made of a Re/Au bilayer (with transition temperature T c pa 1.1 K) 
in the case of a massive rope, and of a Ta/Au bilayer (with T c ph 0.4 K) in the case of a 
thin rope. The room-temperature resistance of the ropes was consistent in each case with 
a resistance of the constituent metallic nanotubes below h/e 2 ~ 25.8 kO [1.1] , which is of 
the order of the value h/Ae 2 corresponding to ballistic transport in a nanotube. 

Below the transition temperature of the contacts, the nanotubes can support currents 
without developing any resistance [TJ- These so-called supercurrents reach a maximum 
value I c , which is known as the critical current. This has been measured for the samples 
described in Ref. [TT] and, as pointed out there, the magnitudes and temperature depen- 
dences obtained cannot be understood in the framework of the conventional proximity 
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effect. The critical current is usually related to the values of the normal resistance Rn 
and the energy gap A in the superconductor through the expression J c = (ii/2)A/eRN- 
However, the value estimated in this way turns out to be about 40 times smaller than the 
experimental measure for the thin rope reported in Ref . ^T] . The observed dependence of 
the critical current on temperature is also very flat in that sample, before getting close to 
the transition temperature of the electrodes. This is again at odds with the conventional 
picture, in which the critical current should reflect the temperature dependence of the 
superconducting gap at the contacts. 

A consistent explanation of the behavior of the supercurrents measured experimentally 
has been proposed in Ref. [TH]. It has been stressed there the need to consider the 
supercurrents as an effect of one- dimensional (ID) transport along the carbon nanotubes. 
The propagation of the Cooper pairs in the nanotubes is favored in the thick ropes, as there 
is a larger attenuation in the strength of the Coulomb interaction for higher content of 
metallic nanotubes in a rope The origin of this effect is similar to that of screening in 
a three-dimensional (3D) conductor. In a rope, however, the tunneling amplitude between 
neighboring metallic nanotubes is highly suppressed in general, due to the misalignment 
of the respective carbon lattices [TTl I18j. The absence of a significant intertube electron 
hopping is what keeps the Coulomb interaction long-ranged in the nanotubes, despite the 
large reduction of its strength in the thick ropes. 

In the present paper, we confront the results from a microscopic model of the nanotube 
ropes with the experimental measures from Ref. ^1] . The propagation of the Cooper pairs 
along the carbon nanotubes gives rise to supercurrents with a strong dependence on the 
length of the junction, when the Coulomb interaction prevails in the electron system. 
Thus, the measures of the supercurrents, available in two experimental samples with 
different length and content of metallic nanotubes, can be used to discern the character of 
the interactions present in the nanotubes. We will see, in particular, that the value of the 
critical current in the thin rope reported in Ref. JT] can be accounted for by considering 
an attractive component in the electron-electron interaction, coming from the coupling of 
the electrons to the elastic modes of the carbon lattice. In general, the balance between 
the renormalized Coulomb interaction in the ropes and the attractive component of the 
interaction allows to explain the magnitudes as well as the temperature dependence of 
the critical currents reported in Ref. [TTj . 

In the next section, we review the different sources of electron-electron interaction 
in the carbon nanotubes and describe the microscopic model relevant for the nanotube 
ropes. Section 3 is devoted to set up the framework used for the study of transport along 
the junctions. The comparison between the results from the microscopic model and the 
experimental data is carried out in Section 4. Finally, the conclusions of our investigation 
are drawn in Section 5. 

2 Electronic interactions in carbon nanotubes 

We review first the structure of the modes and interactions that take place at low energies 
in the carbon nanotubes. The band structure of all the metallic nanotubes has two pairs 
of linear branches that cross at opposite Fermi points (for undoped nanotubes) as shown 
in Fig. [TJ For typical nanotubes in a rope with radius R around ~ 0.7 A, the electron 
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Figure 1: Left: Schematic representation of the wrapping action leading to an armchair 
nanotube. Right: Band structure of an armchair nanotube with nine subbands. The 
energy is measured in units of the overlap integral and the momentum in units of the 
inverse of the C-C distance. 

modes in the energy range below « 0.5 eV belong to the linear branches near the Fermi 
level, as can be seen from Fig. ^ Therefore, we will pay attention only to them when 
studying the low-energy properties of the carbon nanotubes. 

We will carry out the discussion throughout this section having in mind the specific 
geometry of the armchair nanotubes, illustrated in Fig. ^ while a completely parallel 
construction can be done for the metallic chiral nanotubes. The linear branches found 
near the Fermi level can be divided in two classes, according to the different symmetry of 
their modes under the exchange of the two sublattices of the hexagonal carbon lattice. In 
the case of armchair nanotubes, the two subbands that cross at the Fermi level, labeled 
1 and 2 in Fig. Q correspond to the eigenvalues of the one-particle hamiltonian 

/ 1 - 2cos( v / 3A;a/2) 

^ 1 -2cos( v / 3A;a/2) 

that operates in the space of the electron amplitudes in the two different atoms of the 
lattice basis[19j. The parameter t represents the matrix element of the potential for 
nearest-neighbor carbon orbitals and a is the C-C distance. Then, it is clear that the 
modes in subband 1 have a smooth amplitude in the carbon lattice, while the modes 
in subband 2 have an amplitude that alternates the sign when shifting between nearest- 
neighbor atoms. This has important consequences regarding the form of the interactions 
between the different branches. 

2.1 Coulomb interactions 

The Coulomb interaction provides a strong source of repulsion between the electrons 
in the carbon nanotubes. One may obtain an effective ID potential Vc{x — x') on the 
longitudinal dimension of the nanotube, by taking the average of the 3D interaction over 
two nanotube sections placed at a spatial separation x — x' [H]. The Fourier transform 
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with respect to this variable gives a potential which is logarithmically singular for small 
momentum transfer k [20J, 

V c {k) = {e 2 /2n)log\(k c + k)/k\ (2) 

The large momentum k c , of the order of ~ 1/-R, survives in the projection onto the 
longitudinal dimension as a remnant of the finite radius of the nanotube. 

Actually, the Coulomb interaction remains singular only when acting between currents 
that belong each to a given linear branch. Otherwise, if one of the electron modes is 
scattered to another branch with different symmetry in the two-atom basis, the product 
of the in and out amplitudes has an alternating sign over the hexagonal carbon lattice. 
Then, the scattering amplitude averages to zero when making the ID projection with a 
large spatial separation x — x' along the nanotube. 

At short distances, the Coulomb interaction may still give rise to a residual ID interac- 
tion mixing in and out modes in different linear branches, but these processes do not bear 
the singular character of the long-range interaction. This also applies to the processes 
with large momentum transfer in which the electron modes are scattered from one of the 
Fermi points to the other. These backscattering and Umklapp processes have a strength 
which is reduced by a factor of the order of ~ 0.1a/ R with respect to the nominal strength 
e 2 /27r of the long-range Coulomb interaction [Sj. For nanotubes with R « 0.7 A, that 
factor has a magnitude of ~ 0.01. Thus, the backscattering and Umklapp interactions 
of this kind turn out to be negligible in comparison to the Coulomb interaction between 
currents in definite linear branches. 



2.2 Electron-phonon interactions 

Another source of electron-electron interaction comes from the coupling of the electrons to 
the elastic modes of the carbon lattice. The exchange of phonons gives rise to a retarded 
interaction, which can be represented in terms of an effective potential V(k, u) of the form 

V(k,u) = -9 P ,p'(k)g q , q >(-k) — (3) 

— UJ -f- UJ k 

where uj^ is the phonon energy and the g PtP '(k) are appropriate electron-phonon couplings. 
Thus, for frequencies u below the characteristic phonon energies, the potential V(k,u) 
may reach negative values, providing then a kind of effective attraction between the 
electrons. 

The strength of the effective interaction depends on the electron-phonon couplings 
9p,p'{k), where the indices p and p' label the gapless subbands to which the in and out 
electron modes belong. For the computation of the g p , p '(k), we rely on the tight-binding 
approximation, which is fairly appropriate for the carbon nanotubes. The electron-phonon 
couplings can be represented by a sum over nearest neighbors of the atoms in the unit 
cell of the nanotube [21 

g PiP ,(k - k') = l —- Tf2 J2 u^{k)u { p{k'){e s {k - k') - e s ,{k - k')).VJ(s, s') (4) 

where e s (k — k') is the phonon polarization vector at site s, u^\k') and u^(k) are the 
respective amplitudes of the incoming and outgoing electrons, J(s, s') is the matix element 
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of the potential between orbitals at s and s', and fj, is the mass per unit length. We are 
already implying that the coupling can be expressed as a function of k — k', what is a 
reasonable approximation for small values of that momentum transfer (modulo 2kp). 

The product of the two couplings in (jHJ) can be actually positive or negative, depending 
on the linear branches involved[21J. The sign is dictated by symmetry rules which are 
direct consequence of the geometry of the nanotube lattice. We illustrate that fact for the 
case of the armchair nanotubes, though a similar argument can be developed in general 
for metallic chiral nanotubes. Thus, taking into account that the modes with antibonding 
character have an amplitude with alternating sign in nearest-neighbor sites, we observe 
from inspection of (0J) that 

9i,i(k) = -g 2 ,2(k) (5) 

The electrons can be also scattered from one of the subbands to the other by exchanging 
phonons. The fact that one of the modes has then antibonding character reflects in the 
antisymmetric condition 

9i,2(k) = -92,i(k) (6) 

The relations (0) and (jHJ) are based on symmetry rules, and they hold irrespective of 
whether the exchanged phonons are optical or acoustic. 

Using the fact that g* p i{k — k') = g p >, p (k' — k), it is clear from (J3J) that the exchange 
of phonons gives rise to an attractive electron-electron interaction when the electrons are 
scattered within the same subband. On the other hand, the condition (jSJ) on the couplings 
makes the interaction repulsive between electron currents in different subbands. These 
considerations hold for frequencies below the typical Debye energy of the phonons, which is 
in the range between 0.1 and 0.2 eV for nanotubes with radius around ~ 0.7 A[22[I23I21]- 
Taking into account these magnitudes, one can make an estimate of the strength of the 
effective interaction, which turns out to be of the order of ~ [dtj dd) 2 j ~ O.Ivf- 

2.3 Low-energy dynamics 

When describing the low-energy dynamics of the electrons, one has to bear in mind the 
special structure of the nanotube ropes. These are made of a disordered mixture of 
nanotubes with different helicities and diameters. In these conditions, the tunneling of 
electrons between neighboring metallic nanotubes is highly suppressed. This comes from 
the fact that the carbon lattices of neighboring nanotubes are not aligned in general, what 
leads to a mismatch between the respective Fermi points and to the consequent difficulty 
to conserve the longitudinal momentum in the tunneling process[T7]. An experimental 
signature of the very small tunneling amplitude has been observed in the measurement of 
the coupling resistance between tubes, which has shown wide variations ranging between 
2 MQ and 140 MQ in different samples [TH]. This is consistent with a picture in which the 
tunneling amplitude between nanotubes may be about three orders of magnitude below 
the electron hopping amplitude within the carbon lattice. We will therefore assume that 
the coupling between nanotubes is given by the Coulomb interaction, and that the kinetic 
energy of the electrons can be represented by the sum of the kinetic energies of the metallic 
nanotubes in the rope. 

Within each nanotube, one has to balance the repulsive Coulomb interaction and 
the effective electron-electron interaction coming from phonon exchange. As remarked 
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Figure 2: Linear branches and respective electron fields in the low-energy spectrum of a 
metallic nanotube. 

above, there are no significant backscattering or Umklapp processes coming from the 
Coulomb interaction, so that the competition between that and the effective attraction 
takes place in processes where the scattered electrons remain in their respective branches. 
Our approach will focus on the description of this balance between repulsive and attractive 
interactions. For that purpose, we will build a model in which the elementary objects are 
the electron density operators, given in terms of the Fermi fields W^ a {x) for the different 
linear branches (shown in Fig. EJ) by 



The label a runs over the different metallic nanotubes in a rope, a = 1, . . .n. The index 
r = ± is used to label the left- or right-moving character of the linear branch, and the 
index % = ± to label the Fermi point. In this way, we are setting aside the backscattering 
processes, which tend anyhow to enhance the superconducting correlations |25j. Then, it 
can be thought that our model leads to a slight underestimation of these effects, for a 
given strength of the bare attractive interaction. This can be otherwise corrected by a 
suitable renormalization of the effective coupling derived from phonon exchange, assuming 
that it reaches greater values in the low-energy regime of the model[26j. 

According to the above considerations, we write the hamiltonian for a nanotube rope 
in the form 



The first term in (JHJ) comes from the known representation of the kinetic energy for a 
ID electron system [2*7} I2~H] . The potential V^J(fc) mediates the interactions within a 
given nanotube for a = b and between different nanotubes for a / 6. Due to the long- 
range character of the Coulomb interaction, its contribution has to be included in all the 
elements . It is only for a = b, however, that one has to add the contributions arising 
from phonon exchange within each nanotube. 

We can introduce two plausible assumptions that lead to a straightforward resolution 
of the model under study. The first is that all the interactions in the nanotube system are 
spin-independent. Then, it is appropriate to introduce the symmetric and antisymmetric 
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combinations of density operators for the two spin projections. Only the charge density 
operators 

plt(k) = ^(£\(k)+p%(k)) (9) 

enter in the interaction term of the hamiltonian (|Sj) . This implies, in particular for the 
calculation of the correlators, that half of the degrees of freedom in the model are not 
affected by the interaction. 

We will further assume that the strength of the effective interaction from phonon 
exchange does not depend on the particular linear branches to which the scattered electron 
modes belong. We know however that, if the interaction takes place between currents in 
different subbands, the condition (|HJ) introduces a relative minus sign. It is convenient 
to define the respective charge density operators for the bonding and the antibonding 
subband 

= $U k ) + p-U k ) ( 10 ) 

Pt{k) = p^l p {k)+p%{k) (11) 

The Coulomb interaction has always repulsive character, and it acts on the symmetric 
combination of (|10p and (fTTj) 
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The effective interaction from phonon exchange is attractive between currents in the 
same subband, but repulsive for currents in different subbands. It therefore acts on the 
antisymmetric combination 

$(k) = ±($(k)-$(k)) (13) 

In the new density variables, the hamiltonian can be written in the form 

1 r k 



Ho = 2 Vp L dk Y, '■ Pril(k)p { r i(-k) : 

' Cc aria 

+\ [ kc ^ ( A T,p { :l(k) vfc(*)£ ~p+i(-k) 

z J~k c i-k \ a b 

+AgY,P^(k)p { -l(-k)) (14) 

a ) 

where g (< 0) parametrizes the strength of the effective interaction from phonon exchange. 
The hamiltonian (fT^j) encompasses the main interactions in the nanotube system, provid- 
ing a suitable starting point to study the superconducting correlations in the nanotube 
ropes. 



3 Superconducting correlations and supercurrents in 
carbon nanotubes 

We introduce in what follows the theoretical framework appropriate for the description of 
supercurrents in the carbon nanotubes. The study of the Josephson effect in a Luttinger 
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liquid placed between macroscopic superconductors has been carried out in Refs. and 
[3*U] . We first review briefly the approach proposed in Ref. [3U] for the calculation of the 
critical current through that kind of junction. 

In order to describe the coupling to superconducting electrodes, one more term H T j 
has to be added to the hamiltonian of the ID system for each of the tunnel junctions 
created at the contacts. The term Hxj expresses the hopping of electrons from the ID 
system to the superconductor and viceversa, with a given tunneling amplitude tj. It has 
then the form 

H Tj = tj £ ^(d^sjAdj) + h.c. (15) 

where ^/sj,a(dj) is the electron field at the end of the superconductor Sj, and ^/^\dj) is 
a generic electron field of the ID system. 

On top of that, one has also to include the hamiltonians describing the condensates 
for the respective macroscopic superconductors. For our purposes, all the relevant infor- 
mation about their properties can be given in terms of the respective order parameters 
Aj and normal densities of states Nj at the contacts. 

The supercurrent I through the ID system can be computed in terms of the derivative 
of the free energy T with respect to the phase difference \ °f t ne order parameters [30 

dT 

I( X ) = -2e— (16) 
OX 

When the junctions are not highly transparent, one can adopt a perturbative approach 
by expanding in powers of the tunneling hamiltonians Htj- For usual instances in which 
the time of propagation between the tunnel junctions is larger than 1/A, the contribu- 
tions to the derivative in (|T6*j) are dominated by processes in which the Cooper pairs 
tunnel from one of the superconductors to the ID system, propagating then to the other 
superconducting contact [SUl- These first perturbative contributions give the result 

I(X) ~ 8etlt 2 2 k B T^- f 1/kBT dn r dr 2 H dr 3 H dr±(¥ si a (d u -in)^ si _ a (d u -ir 2 )) 
d\ Jo Jo Jo Jo 

<*W(d!, -in)*to( dlt - iT2 )^)\d 2 , -ir 3 )¥$(d 2 , -<r 4 )> 

(*S2, CT '(d 2 , -iT 3 )y S 2,-*i(d 2 , -ir 4 )) (17) 

where the statistical averages, at temperature T, are taken over ordered products with 
respect to imaginary time r, in the respective systems of the fields. 

At low enough temperatures, we can approximate the statistical averages at the bound- 
ary of the superconductors by delta functions in imaginary time 

(*s J> (d i , -in^sjMdj, -ir 2 )> ~ e^A^n - r 2 ) (18) 

Xj being the phase of the order parameter. Then the expression (fTTjl is greatly simplified. 
The maximum value of the supercurrent, that we call critical current I C {T), is given by 

I c {T)^2eN x N 2 t\tl /^V^d!,-^)^ (19) 
Jo 

We now make the passage to apply the above approach to the case of the carbon 
nanotubes. The critical current in ()19j) is given in terms of the propagator of the Cooper 
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pairs along the ID system. We need therefore to compute that object with the fields 
corresponding to the linear branches crossing the Fermi level of the carbon nanotubes, in 
the model governed by the hamiltonian (|14j). 

The calculation of the ID correlators can be accomplished by making use of bosoniza- 
tion techniques |27| For each linear branch of the carbon nanotubes, one can define a 
boson operator related to the respective electron density operator by 

0,*&(aO = 2*p<&(x) (20) 

The bosonization method relies on the feasibility of expressing the fermion fields in terms 
of the respective boson operators: 

= A= ex P Mit(z)) (21) 

where a c is a short- distance cutoff of the order of k^ 1 |27j . 

In the case of the carbon nanotubes, the propagator that enters in the computation 
of the critical current in (|19|) is of the type 

G{x,t) = (^V ' )*-;^ )*-;^*)*^^.*)) ( 22 ) 

where the choice of linear branches accounts for the fact that the Cooper pairs are formed 
with zero total momentum. The evaluation of (|22|) can be done by using the bosonization 
formulas (|2*T|) and passing to the combination of boson fields that bring the hamiltonian 
(fTH) into diagonal form. 

We observe then that, in the nanotube ropes, the Coulomb interaction acts only on the 
combination of density operators corresponding to the total charge density in the rope. 
That is, the repulsive interaction becomes sensible only in one out of the An possible 
interaction channels, arising from the n different metallic nanotubes and the degeneracy 
in the gapless subbands and spin projections. 

The An different channels can be classified into the sectors corresponding to the sum 
and the difference of electronic charge in the gapless subbands, and the rest which take 
into account the spin densities for the two spin projections. The expression of correlators 
like (}2"2"|) factorizes into contributions for each of the different sectors, in such a way that 
the propagator for the Cooper pairs can be written in the form 

i n 3n— 1 

G(x, t) = —C(x, t) J] N(x, t) J] F(x, t) (23) 

a c 1 1 

The first factor comes from the propagation of the total charge density from all the 
metallic nanotubes, which has a purely electrostatic interaction. The next n equivalent 
factors stand for the contribution of the charge mismatch in the gapless subbands, which 
bears the effect of the attractive interaction according to (|14|) . The remaining 3n — 1 
factors account for the free propagation of the rest of density excitations, including the 
spin degrees of freedom and charge degrees of freedom orthogonal to those affected by the 
interaction. 
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Each of the factors in ()23|) corresponds to the propagation of a boson field. At zero 
temperature, for instance, all of them adhere to the common expression 

X(x,t) = expi — — [ C dk ,\ , (1 - cos(kx) cos(v F kt))\ (24) 
\ 2n Jo n{k) k J 

The respective quantities n(k) arise in the change to the variables that make the hamil- 
tonian (JT4*j) diagonal. This process also leads to a renormalization of the respective 
Fermi velocities, given by v F (k) = v F /fi(k). In the case of C(x,t), we have //(&) = 
l/\Jl + 4nVc(k)/iiVF, while for N(x,t) the quantity fi(k) does not depend on the mo- 
mentum, ji = 1/ \Jl — 4\g\/nv F . For the rest of excitations with free propagation, we 
simply have p,(k) = 1 and v F (k) = v F . 

In practice, one is interested in the behavior of the supercurrents at finite temperature, 
which requires moreover the knowledge of the propagator (}2*2*j) at imaginary time. For this 
reason, it is more appropriate to deal with the representation of objects like (}2"4*|) at finite 
temperature in the Matsubara formalism, introducing the sum over discrete frequencies 
uj m = 2nmkBT : 

(I fk* 2k B T "^+°° 1 - cos(A:x) cos(u m r) \ 

X(x, -it) = exp -— / dk ^ 25 

V 2n Jo v F m ±^ 0O (uj m /v F y + k 2 J 

Thus, the framework developed is well-suited to discuss the behavior of the supercur- 
rents in carbon nanotubes depending on the temperature T, the distance d\ — g?2 between 
the superconducting contacts, the number n of metallic nanotubes in a rope, and the 
competition between the repulsion given by the Coulomb interaction and the effective 
attraction given by phonon exchange. 



4 Comparison with experimental data 

In this section we establish the comparison with the experimental measures of supercur- 
rents in carbon nanotubes reported in Ref. ^T]. The supercurrents were observed in two 
nanotube samples with different structure. One of them consisted of a rope containing 
approximately 200 nanotubes, with a length of about 1.7 /im. The other sample was made 
of seven tubes merging at one of the ends into a single nanotube, extending to a total 
length of 0.3 fim. It has been remarked in Ref. ^T] that the critical current I c ~ 0.1 fiA 
measured in the thin rope was about 40 times larger than expected within the conven- 
tional framework of the proximity effect. The very flat shape of the critical current as 
a function of the temperature was also surprising, as its decay was only apparent in the 
neighborhood of the critical temperature for the superconducting electrodes. The thick 
rope showed a more conventional behavior, although the decay of the supercurrent at the 
critical temperature of the electrodes was not as sharp as expected [TT]. 

We face the challenge of accounting for the experimental values and features of the 
critical currents with our microscopic model, that allows to describe ropes with any kind of 
nanotube content. In this respect, our computational framework has only a few adjustable 
parameters. The prefactors that affect the strength of the critical current in (|19|) can be 
encoded into the relative conductances Gi at the tunnel junctions, given (in units of 4e 2 ) 
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Figure 3: Plots of the critical current (measured in units of ev F k c ) versus distance between 
the superconducting contacts (in units of k" 1 ) at T = 0. The full lines represent the results 
for ropes with a number of metallic nanotubes n = 60 (top) and n = 2 (bottom), taking 
2e 2 /ir 2 v F = 1.0 and a strength of the effective attractive interaction g/nv F = 0.2. The 
dashed lines represent the respective counterparts for the same values of n in the absence 
of attractive interaction. 



by Gi = Nit 2 /vp- The critical current gets the right dimensions from the combination 
ev F k c , that arises after setting a" 1 = k c and performing the integral in (|T§|) over the 
dimensionless variable f = v F k c r. The expression of the critical current becomes 

r v F k c /k B T n 3n-l 

J C (T) « 2G 1 G 2 ev F k c / d?C(L, -if/v F k c ) J[ N(L, -if/v F k c ) f[ F(L, -if/v F k c ) 

Jo 11 

(26) 

where L represents the distance between the superconducting contacts. 

Equation (J26)) gives the critical current per metallic nanotube in a rope. We observe 
that the order of magnitude is given at short distances by the quantity ev F k c . It is 
reasonable to set the short distance cutoff according to the scale of the nanotube diameter, 
and we will take k c « 0.25 nm" 1 for typical nanotubes in a rope. This sets a unit scale 
ev F k c »s 15 fiA for the critical currents in carbon nanotubes. 

One has to bear in mind, however, that the critical currents obtained from (|26|) display 
a strong dependence on the distance between superconducting contacts. Their decay for 
increasing length L can be very pronounced in ropes with a small number of metallic 
nanotubes, where the electron interactions are predominantly repulsive. The behavior 
at large L becomes much softer in the thick ropes, as shown in Fig. El In general, the 
decay of the critical currents with the length L is not given by a simple power-law, since 
the dependence of the Coulomb potential Vc(k) on the momentum prevents the perfect 
scaling of the observables at long distances. 

The results shown in Fig. |3] allow us to draw the comparison with the values of the 
critical currents at T rs K reported in Ref. [llj. The appropriate correspondence 
with the thin rope described there is established by making the reasonable assumption 
that two out of the seven nanotubes in the rope are metallic. Moreover, we determine the 
theoretical value of the critical current by looking at the decay along the distance between 
the contacts, setting L = 75/ k c for the thin rope. For that choice of the parameters 
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and taking a strength 2e 2 /ir 2 VF = 1.0 for the Coulomb interaction, we observe that the 
experimental value of the critical current (I c « 0.1 //A) can be matched by taking a 
coupling for the attractive interaction g/itVF = 0.2. Actually, the value of the integral in 
()26|) at zero temperature for n = 2 and L = 75/ k c is ~ 3.8 x 10~ 3 , which gives a critical 
current 2/ c (0) ~ 0.11 //A (times 2G1G2) after adding the contributions of the two metallic 
nanotubes. 

The precise determination of the critical current is affected by the values of the relative 
conductances G\ and G<i- It is however reassuring the fact that, once the interaction 
strengths are fixed to the above mentioned values, one can reproduce approximately 
the ratio between the critical currents observed in the two different samples where the 
supercurrents have been measured. In the case of the thick rope, we assume that it may 
have a number of metallic nanotubes n ~ 60, and we make the estimate of the critical 
current by computing the decay for a 1.7 /im-long rope, which corresponds to L — 425/ k c 
in our model. The value of the integral in ()26|) at zero temperature for n = 60 and 
L = 425/ k c is ~ 5.1 x 10~ 3 . This leads us to assure that the critical current in the thick 
rope should be about 40 times greater than that in the thin rope, taking into account the 
approximate ratio between the number of metallic nanotubes in the two samples. This 
agrees well with the order of magnitude of the experimental value I c ~ 2.5 fiA for the 
thick rope reported in Ref. jllj . 

The other important check is that our model is able to account for the dependence of 
the critical currents on temperature observed experimentally. We have represented in Fig. 
El the behavior obtained from (|26p. computing now the critical current at a fixed distance 
L = 75/ k c , for the thin rope with n = 2 (and keeping a suitable value of the coupling 
g/nvp = 0.2). The temperature is measured in units of vpk c /ks ~ 1.2 x 10 3 K. Thus, in 
the case of the thin rope, the transition temperature T c m 0.4 K of the superconducting 
contacts corresponds to a value k B T c /vpk c ~3x 10~ 4 in the scale of Fig. |U We observe 
that the natural behavior of the critical current is being flat over such range of small 
temperatures, in agreement with the experimental measures in the thin rope studied in 
Ref. [TT]. The present results can be trusted while remaining away from the critical 
temperature T c of the superconducting electrodes, as we have not paid attention to the 
dependence of their gaps on temperature. The sharp decay of the critical current observed 
experimentally near T c is just the natural consequence of approaching the superconducting 
transition of the contacts. 

We have to remark, however, that the flat behavior of the critical current as a function 
of T is a consequence of the short distance along which it has been measured in the thin 
rope. For long enough ropes, there is a crossover temperature above which the behavior 
changes to a more pronounced decrease of the critical current. This can be appreciated 
in Fig. 01 where it has been represented the critical current that would correspond to 
a thin rope with n = 2 and length L = 425 /k c . In that case, the critical current at 
zero temperature would be 2J C (0) ~ 0.009 fiA (times 2G\G2)- Most remarkably, an 
inflection point is observed in the dependence on temperature. On physical grounds, such 
a crossover temperature marks the point above which the thermal effects begin to affect 
significantly the propagation of the Cooper pairs, inducing the pronounced decrease of 
the super current. 

The crossover temperature for the critical current is always found below a scale which 
is of the order of ~ vp/ksL. This is perfectly consistent with the behavior of the 1.7 ^m- 
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Figure 4: Plots of the normalized values of the critical current versus temperature (mea- 
sured in units of vf^cI^-b) in a thin rope with 2 metallic nanotubes, for 2e 2 /% 2 vf = 1.0, 
g/nvF = 0.2, and respective nanotube lengths L = 75/ k c (full line) and L = 425/ k c 
(dashed line). 

long rope described in Ref. [H], where an inflection point can be seen in the plot of 
the critical current as a function of the temperature. We have represented in Fig. 03 the 
results for the critical current computed from (|26jh for n = 60 and L = 425 /k c . Bearing 
in mind that the temperature is measured in the figure in units of vphc/kg ~ 1.2 x 10 3 K, 
we observe that the inflection point in the curve of I C {T) corresponds to a crossover 
temperature T s; 0.7 K. This is in close agreement with the experimental measures 
reported in Ref. [TT] . 

We have also represented in Fig. El the behavior that should be expected if the thick 
rope were shrinked down to a length of ~ 0.3 /im, keeping the same number of metal- 
lic nanotubes. In that case, the critical current at zero temperature could be as large 
as n/ c (0) ~ 22.5 /xA (times 2GiG 2 ). In agreement with our preceding discussion, the 
crossover temperature would be shifted in the shorter sample to T ~ 5 K, a value that 
would be in general larger than the transition temperature of the contacts and would 
make then the crossover unobservable. 

5 Conclusions 

In this paper, we have developed a transport approach with the aim of accounting for the 
properties of the supercurrents measured in carbon nanotubes. Some important features 
of the critical currents reported in Ref. ^T] cannot be explained within the conventional 
picture of the proximity effect. In that framework, the values of the supercurrent at zero 
temperature are given by the expression irA/eR^, where Rn is the normal resistance 
of the junction. That estimate falls short by a factor of 40, for instance, to reproduce 
the critical current measured in the 0.3 /im-long rope considered in Ref. [11J. In the 
ID approach we have adopted, the values of the critical currents are determined instead 
by the scale of the distance between the superconducting contacts, provided that such 
a quantity is greater than the coherence length in the superconductors. The theoretical 
values found in this way are in agreement with the experimental measures for the two 



13 



0.2 



0.00025 0.0005 0.00075 0.001 
T 

Figure 5: Plots of the normalized values of the critical current versus temperature (mea- 
sured in units of df^c/^b) in a thick rope with 60 metallic nanotubes, for 2e 2 /h 2 vf = 1.0, 
g/itVF = 0.125, and respective nanotube lengths L = 425/ k c (full line) and L = 75/ k c 
(dashed line). 

samples where supercurrents have been reported in Ref. |llj . 

Even in a ID framework, it is possible to make two different descriptions depending 
on whether the contacts are supposed to be perfectly transmitting or not. In the first 
case, the decay of the critical current turns out to be independent of the strength of the 
interactions at low temperatures |2Sj. Such an instance may not be the most appropriate to 
describe the experimental conditions of the measurements of the proximity effect described 
in Refs. [TT] and ^2]- In the latter, no evidence of supercurrents has been found for 
individual nanotubes with a length of about 0.3 /im, placed between Nb contacts. One 
may think that the main difference between the conditions of the experiments described 
in Refs. []' and lies in that, in the former, the samples were made of ropes suspended 
between the electrodes. Such conditions seem to be determinant for the observation of 
supercurrents, what is supported by the fact that the normal resistances of the junctions 
were comparable in the two kind of experiments. This should lead us to conclude that the 
electron-electron interactions play an important role in the observation of the proximity 
effect in the mentioned experiments. 

We have consequently adopted the point of view in which the single-particle scattering 
dominates at the contacts between the nanotubes and the electrodes. In these conditions, 
the values of the critical current depend sensibly on the strength of the Coulomb repulsion 
and the attractive interaction arising from phonon exchange. In the thinner ropes, the 
repulsive interaction does not suffer a significant reduction from the electrostatic coupling 
among a small number of metallic nanotubes, and this explains the comparatively small 
magnitudes of the supercurrents in those instances. The value of the critical current J c w 
0.1 fiA for the very thin rope reported in Ref. ^T] is actually dictated by the relatively 
small distance (~ 0.3 /im) between the contacts in that sample. We can foresee that, in 
a similar rope with only two metallic nanotubes extending over a distance L ~ 1.7 ^m 
between superconducting contacts, the value of the critical current should decay to a 
magnitude of ~ 10~ 2 /iA. 

In the thicker ropes, larger supercurrents are expected as the strength of the Coulomb 
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repulsion is significantly reduced from the interaction among the different metallic nan- 
otubes. On the other hand, it is known that the coupling resistance between tubes in a 
rope is typically above 1 Mfl, what leads to assume that the conduction is confined to 
the individual nanotubes. Then, each metallic nanotube in the rope provides a different 
channel for the propagation of Cooper pairs, whenever there is good contact with the su- 
perconducting electrodes. This condition can be difficult to ensure from an experimental 
point of view, introducing some uncertainty about the real number of metallic nanotubes 
that may contribute to carry the supercurrent. We have seen, in particular, that it is 
possible to account for the value of the critical current in the thick rope described in Ref. 
[TT] . by adding the currents for the approximate number of metallic nanotubes in the 
rope. 

Finally, let us remark that we have also given a satisfactory description of the temper- 
ature dependence of the critical currents reported in Ref. jllj . The steady behavior of the 
critical current in the thin rope studied there can be understood as a consequence of the 
relatively short distance between superconducting contacts for that sample. We have seen 
that, for sufficiently long nanotubes, it is possible to observe a crossover in temperature 
from the flat behavior to a pronounced decay of the critical current, before the transition 
temperature of the contacts is reached. The crossover has to take place at a temperature 
which is in general slightly below the scale of vp/ksL, and it manifest by the presence 
of an inflection point in the plot of the critical current as a function of the temperature. 
Such a feature has been actually seen in the experimental measures of the 1.7 /im-long 
rope studied in Ref. jTTj. This comes as a genuine signature of the propagation of the 
Cooper pairs along the nanotubes, and stresses once again the suitability of the transport 
approach to account for the unconventional properties observed in the carbon nanotubes. 
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